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1. Introduction 


Scattering amplitudes of gluons in Yang-Mills theories exhibit a remarkable simplicity 
that is not manifest from their calculation using Feynman diagrams. At tree-level, Parke- 
Taylor or maximally helicity violating (MHV) amplitudes [Q provide a striking example. 

One-loop MHV amplitudes in A/" = 4 super Yang-Mills also exhibit remarkable sim¬ 
plicity when computed using the unitarity based method . Recently, a new technique 

was presented for computing general one-loop amplitudes in A/” = 4 super Yang-Mills using 
quadruple cuts p. This is a systematic and simple procedure that blends old |P,|^, more 
modern B»B0ra and very recent ideas to uncover the simplicity of generic 

amplitudes. Using this one can easily reproduce all known results PJ^, pR| , p^ , pR| , p!^ and in 
principle compute any other amplitude. 

One main motivation for computing A/” = 4 amplitudes of gluons is that they are part 
of amplitudes in theories with less supersymmetry |P] (for a review, see [^|). 


In this paper we concentrate on one-loop Af = 1 amplitudes of gluons. Such an 
amplitude can be decomposed as follows. 




J\f=l vector 




A/'=4 


- 3A 


Ar=i 


( 1 . 1 ) 


where is an amplitude where the full J\f = 4 multiplet runs in the loop, and A^^^ 

denotes the contribution from an A/” = 1 chiral supermultiplet running in the loop. As 
mentioned above, the Af = 4 problem is easy to solve using quadruple cuts. On the other 
hand, A^^^ only contains fermions and scalars in the loop and thus it is expected to be 
simpler than the full Af = 1 vector multiplet. This is why the decomposition (|1 . Ij ) is useful. 

The computation of A^^^ is also important because it is part of a supersymmetry 
decomposition of a QCD amplitude at next-to-leading order. 


^QCD _ 


A/'=4 


4A 


M=l 




scalar 


( 1 . 2 ) 


where denotes an amplitude with only a gluon running in the loop. is an 

amplitude with only a complex scalar running in the loop. 

The beneht of this approach is that supersymmetric amplitudes are four-dimensional 
cut-constructible [|],^. This means that they can be completely determined by studying 
their hnite unitarity cuts, and therefore the dimensional regularization parameter can be 
set to zero. Furthermore, they can be expressed as a linear combination of known integrals 
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called scalar box, triangle, and bubble integrals, with rational coefficients in the kinematical 
invariants. 

The scalar part is more complicated, as only part of it can be determined by studying 
four-dimensional unitarity cuts. There are single-valued pieces that have to be determined 
using some other method.0 The current state of the art in QCD is the five-gluon ampli¬ 
tude |T^. This means that even the scalar part has been fully computed for all helicity 
configurations. 

For special helicity configurations much more is known. is known for all 

MHV amplitudes |Q]. Also in |^, the cut constructible part of was given 

for MHV amplitudes where the gluons of negative helicity are adjacent. More re¬ 
cently, the non-adjacent case was computed in [j^ using the techniques of []Tl| , [2^ , pT . 
Also recently, 2“, 3“, 4+, S"*", b"*") was presented in and then extended to 

^A/'=i(i-, 2“, 3“, 4+,..., n+) in [^. Also in [^, the scalar box coefficients of all other 
next-to-MHV six-gluon amplitudes were computed using quadruple cuts 0. 

For apart from the five-gluon case, all amplitudes with at most one negative 

helicity gluon are also known [p4|,^ , |2^ . 


In this paper, we introduce a systematic approach to computing any finite unitarity 
cut of amplitudes of gluons. Although our focus is on it is important to mention 

that this can also be applied to obtain the four-dimensional cut constructible part of 
or even as an alternative way of computing amplitudes. The basic idea is to exploit 

the representation of the Lorentz invariant measure of a null vector £, introduced in ||2CI|| , 
as a measure over R+ x CP^ x CP^ with contour of integration a certain diagonal CP^. 
More explicitly, one writes ^aa = and then 


= / tdt i _(A,dA)[A,dA](«) 


(1.3) 


i\=\ 


where (•) represents a generic integrand. 

It turns out that the integration on the right hand side of (|^) can always be reduced 
in a systematic way to an integral performed in |^. The main simplification arises be¬ 
cause the final integrals always localize to some poles in the region of integration. This 
is reminiscent of the technique developed in [|^ where certain differential operators are 
applied to the cut integral in order to produce a localization via a holomorphic anomaly 


^ In principle the whole scalar part can be computed from unitarity cuts if higher orders in e, 
the dimensional regularization parameter, are kept. (See section 4.4 of [|li 
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T^. Surprisingly, here we find that up to an integration over a single Feynman parameter, 


which is responsible for logarithms, all unitarity cuts localize by themselves without the 
need of a differential operator. 

In the case of unitarity cuts of A/” = 1 amplitudes of gluons, , one expects bubble, 

triangle and box scalar integrals to contribute to a given cut. Remarkably, our procedure 
naturally leads to a clean separation of the three kinds of contributions, allowing for an 
individual calculation of the corresponding coefficients. 

An important simplification in is that one- and two-mass triangle coefficients 

never need to be computed. It turns out that their contributions always cancel against 
singular pieces in box integrals. This leaves us with bubbles, three-mass triangles and 
finite boxes, which we define in detail. Since the coefficient of scalar boxes can easily be 
computed from quadruple cuts one can disregard that piece and concentrate on the 
bubble and three-mass triangle scalar integral coefficients. 

As an application of our technique we compute all A^^^ next-to-MHV six-gluon 
amplitudes. The reason we have undertaken the whole calculation is because six-gluon 
amplitudes in QCD are going to be important for future colliders and our computation 
completes the second piece in ( |1.2|) . It is important to mention that this calculation requires 
the use of tree-level amplitudes of gluons with two fermions or two scalars. Luckily, very 
compact formulas for those amplitudes were derived very recently p8|j2^] by extending the 
techniques of |30| , ^ . 

This paper is organized as follows: In section 2, we define A^^^ in terms of a linear 
combination of scalar integrals. Using its singular behavior we show that only bubble, 
three-mass triangle and finite boxes are necessary. In section 3, we study general unitarity 
cut integrals and present the method for computing them explicitly. For A^^^ we explain 
the way to basically read off the coefficients of bubbles and three-mass triangles. In section 
4, we present the calculation of A^^^ non-MHV six-gluon amplitudes. We also present 
our results for the next-to-MHV n-gluon amplitude where all negative-helicity gluons are 
consecutive. This amplitude has appeared in ||2^, but by our procedure it emerges in a 
different form. Explicit details are given in order to illustrate the steps described in section 
3. In section 5, we summarize the results for all the amplitudes computed in section 4. 
Section 5 is intended to be self contained so that the reader interested only in the six-gluon 
amplitude results can skip the rest of the paper. Appendix A contains a detailed definition 
of scalar integrals as well as our definition of finite box integrals. Appendix B summarizes 
the results of tree-level amplitudes of gluons and fermions needed in section 4. 
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Throughout the paper, we use the following notation and conventions along with those 
of and the spinor helicity-formalism [p^ ,|33| , ^ . The external gluon labeled by i carries 
momentum Since = 0, it can be written as a bispinor {pi)ad = K aK a- Inner 
product of null vectors Pah = and Qad = can be written as2p - q = (A, A') [A, A'], 


where (A, A') = eaf,A“A'“ and [A, A'] = ea^A^A'”. Other useful dehnitions are: 


Pi...j = Pi +Pi+i H- \-Pj 

= Pi +Pi+1 H- ^Pi+r-l 

tP = {Pi+Pi+1 H-hPi+r-l)^ 

{i\'^Pr\j] = r)[r j] 

r r 

r s r s 

r s r s 

Pt)\j] = r)[r s](s t)[tj] 

r s t r s t 

where addition of indices is always done modulo n. 


(1.4) 


2. One-Loop J\f = 1 Amplitudes 

Amplitudes of gluons at one-loop admit a color decomposition [p5| , |36[| with single and 
double trace contributions. The piece proportional to the single trace term Tr (T“^ ... ) 

is called the leading color partial amplitude and it is denoted by A^.i(l,. .. ,n). In this 
paper we concentrate on A^;i(l,.. ., n). The reason is that when all particles in the 
loop are in the adjoint representation, all sub-leading color amplitudes are given as linear 
combinations of with permutations of the gluon labels (See section 7 of 0 for a proof.) 
This is the case for all amplitudes we consider. In the remainder of the paper we will 
simplify the subscript and just denote the leading color partial amplitude by An{l, ..., n). 

We consider amplitudes of gluons where an A/" = 1 chiral multiplet circulates in the 
loop. Reduction techniques allow us to express these amplitudes in terms of scalar integrals 
in the shapes of boxes I 4 , triangles /a, and bubbles I 2 1^,0- These functions are given 
explicitly in appendix A, along with some helpful hgures. 
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The amplitude thus takes the followiug form: 






(cp/]’" + 4 


2m e 7-2m e i „2m h 7-2m h | ^3m y-Sm i ^4m 7-4m 

A I C/1 -/ /I 


+ crc + crc 


+4""/3^"" + + C 2 / 2 ) 


( 2 . 1 ) 


Here e = (4 — Il)/2 is the dimeusioual regularizatiou parameter, n is the reuormalizatiou 
scale, aud rr is dehued by 

^ rd + c^d-.) 

r(i- 2 c). ' ' 

The sum ruus over all the cyclic permutatious withiu each type of iutegral. The coefhcieuts 
c of the scalar iutegrals are ratioual fuuctious of spiuor products. This follows from the 
reductiou procedure fl. 


2.1. Singular Behavior 

The iufrared aud ultraviolet siugular behavior of these amplitudes is kuowu aud was 
giveu iu p8| , |39| , ^ . For the case of a gluou amplitude with the Af = 1 chiral multiplet iu 
the adjoiut represeutatiou circulatiug iu the loop, the divergeut behavior is giveu simply 
iu terms of the tree-level amplitude by 


■^n I singular 


rr 


-a: 


tree 


where 


rr = 


e(47r)2-^ - , 

r(i + €)rdi-£) 


(2.3) 


(2.4) 


r(i-2£) 

aud is the color-ordered tree-level amplitude. 

Now let us see what this meaus for the scalar iutegral coefhcieuts. The iutegrals /f™' 
aud are huite. Therefore their coefhcieuts do uot coutribute to (|2.3| ). The bubble 
iutegral diverges as 1/e. Oue- aud two-mass triaugle iutegrals cau be couveuieutly writteu 
iu terms of the fuuctiou 


T(,) = J(-s)- 


(2.5) 


as follows: 




1 


T{s), 


1 


(T(s)-T(t)) 


[s) - i-t) 

where s aud t deuote the iuvariauts iu dihereut iudepeudeut chauuels. 


( 2 . 6 ) 
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Finally, one-, two-, and three-mass box scalar integrals have the property that they 
can be made hnite by adding linear combinations of T{s) functions. We denote the hnite 
box integral functions by where F stands for hnite. Their dehnition is given in detail 
in appendix A. 

Now we are ready to derive the main result of this section. From (|2.3D we see that 
all divergences of the form 1/e^ must be absent. This implies that all T{s) functions must 
cancel among the different terms. Since one- and two-mass triangle integrals are given 
entirely as linear combination of T{s) functions with rational coefficients, it follows that 
their coefficients are such that they do not appear in the hnal answer for the amplitude. 
Our interpretation is that the only reason they must be included is to cancel the 1/e^ 
divergences from the box integral. 

Therefore, we reach the conclusion that can be written as a linear combination 

of hnite box scalar integrals / 4 i?, three-mass triangles and bubbles l 2 - More explicitly, 

rvilAy 


A 


N=l 




2-e 


Y, (C2/2 + + c 


Im Tim 


Tim I J2m e T2m e i _ 
J-4F ' ^4 ' ^4 


2m h T2m h 


3m rSm 


L4F 


, 6m T 

I 04 


4F A ^4 

(2.7) 

Finally, among the bubble coefficients there is one relation that must hold in order 
to satisfy namely that the sum of all bubble coefficients reproduces the tree-level 

amplitude: 

= ( 2 . 8 ) 

In section 4, (|2.8|) is used as a very non-trivial consistency check of our results for next-to- 
MHV six gluon amplitudes. 


4m r4m\ 
^4 ) 


3 . Coefficients from Unitarity Cuts 

In this section we introduce a new method for computing explicitly any hnite unitarity 
cut in a gauge theory with massless particles running in the loop. Of course, our aim here 
is to apply the technique to the computation of the coefficients in (|2.7|) which determine 
4 A/'=i Nevertheless, it is important to mention that this can also be applied to obtain 
the four-dimensional cut constructible part of or even as an alternative way of 

computing amplitudes. 

The unitarity cut in the -|- 1,..., j — 1, j)-channel is computed by cutting two 
propagators in the loop whose momenta diher by Pij = pi pj in all Feynman 

diagrams contributing to the amplitude. Adding up all these contribution we hnd a “cut 
integral” 01 


Further information about this technique may be found in Q. This body of work was not 
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Fig. 1: Representation of the cut integral. Left and right tree-level amplitudes are 
on-shell. Internal lines represent the legs coming from the cut propagators. 


f i, i -h 1,..., j - 1, j, + 1, i + 2,..., z - 2, z - 1, 

(3.1) 

where dfi = d^eid^e25^+\il)5^+\il)5{ii + £2 - Pij) is the Lorentz invariant phase space 
measure of two light-like vectors (^ 1 , ^ 2 ) constrained by momentum conservation. See hg. 1. 

This cut integral computes the discontinuity of the amplitude across a given branch 
cut in the space of kinematical invariants. The same discontinuity can be computed from 
the right hand side of O- Note that the coefficients are rational and thus do not have 
branch cuts. The scalar integrals have a cut in this channel only if they contain the same 
two propagators that are cut. 

The idea is to determine the cuts of the known scalar integrals, compute (|3.1[ ) explicitly 
and then solve for the coefficients by comparing both sides. 

The phase space integral in ( |3.1|) can be performed following the techniques of |2^ . 
The idea is to use moment um conservation to write the integral entirely in terms of just 
one of the cut propagator momenta, say £ = ii. This vector is then parametrized as 
£aa = tXa^^cn where the scale t is real and the spinors A and A are independent homogeneous 
coordinates on two copies of CP^. The integral is then performed over the diagonal 
dehned by A = A. The integral can be rewritten as 


(;“«<+>(«^) (•) = / dtt / (A,(iA)lA,(iA](.), 


(3.2) 


intended to apply to massless theories. We find that the material can nevertheless be adapted for 
the considerations of this paper. The modern interpretation is found in 
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where the bullets represent generic arguments. 

We hrst illustrate the procedure by computing the cuts of bubble and three-mass 
triangle scalar integrals entering in (|2.7|) and then we discuss the calculation of the general 
cut integral (p.l|) . 


3.1. Bubble and Three-Mass Triangle Unitarity Cuts 

As a warm-up let us compute the double cut of a bubble and a three-mass triangle 
using (|3.2|) . We will need these results to read off the coefficients we need. Let us denote 
the cut of a scalar integral / by A/. 

The unitarity cut of a bubble is given by (see appendix A for a definition of the original 
integral!) 

Al2{K)= [ - Kf) 


tdt j (A dX)[X dX]6^+\K^ - tKaaX'^X^) 


(3.3) 


= / (A dX)[X dX] 


{KaaX-X-y 

were we have used that in the kinematic regime where > 0 the delta function always 
has its support in the integration region of t. 

We postpone evaluating the last integral in ( ^?^) . First let us evaluate a slightly more 
complicated integral ligi, 


I = {X dX)[X dX] 


(KaaX^^X-) 




with 


g(\) = nfal(A,.4i) 


(3.4) 


(3.5) 


First note the following identity that holds for an arbitrary but hxed negative chirality 
spinor rj: 

v] 


"^^1—^(A) =-dA" ^ 


[KaaX-X-y 


dX^ \{KaaX^X^){KaaX‘^r) 


-d(A) 


(3.6) 


In the following calculations, we are omitting the factor of —i which has been ac¬ 


counted for in the overall factor in ( p.lD and (2.7). 

^ This integral was performed as part of a derivation of MHV diagrams for tree-level amplitudes 


of gluons from a twistor string theory calculation [ H |. 













This identity holds for all values of A except for those where the denominator vanishes 
along the contour of integration. The reason is that along the contour of integration A = A 
and therefore 

-<<vi7^ = 2rf((A,C», (3.7) 

d\^ \A,C) 

where we have introduced a (0, l)-form 5((A, 0), such that 


(A d\) = 


(3.8) 


Let us write the complete form of (|3.6| ) that is valid for all values of A along the 
contour of integration: 


'^^L(7(A) --dA" ^ 


[\v] 


(K„„A“A“)2' 


dX^ \iKaaX-X^)iKaaX<^r) 


-d(A) 


k 




KadX‘^r]^ ^ 
j=i 


(3.9) 

The contribution from the hrst term in (|3.9|) gives zero after integration over A. On the 
other hand, the delta functions in the remaining terms localize the A integral and give the 
value of the integrand at the pole. More explicitly one uses that 


(A,dA)d((A, XB))H{X) = -iH{XB). 


(3.10) 


A short calculation reveals that (|3.4|) is given by 










(3.11) 


where Xk a ^ Ka^rr. 

This is the basic result that will allow us to calculate any double cut in section 3.2. 
Going back to the cut of the bubble integral (|3.3|) we hnd that by setting 5 r(A) = 1 in 

o 

AhiK) = -1. (3.12) 


Consider now a three-mass triangle integral. 
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K, 


K, 


Fig. 2 : A double cut of a three-mass triangle integral. 


Denote the momenta at the vertices by Ki, K 2 and K^. Let ns calculate the cut in 
the Ki channel, Ai. Let the momentum in the cut propagators be £ and i — Ki. See hg. 2. 
Then 


^iir = 


^ ’ {e + Ksr 


'0 J iL|+tiL3,aaA“A“ 


- j {X dX)[X dX] 

- /(A dA)[A dA] 


Kf 


(iLi,„„A“A“ 


(3.13) 


(iLi,„„A“A“)2 iL 32 (iLi,„^A“A“) + KfiKs,aA<^X-) 

1 


(iLi,„^,A“A“)(g„^A“A“) 


where Qad = ^{Ki^ad) + {K 3 ,ad)- 

This integral can be brought to the same form as that in (^) by introducing a 
Feynman parameter to combine the two denominators into one. This Feynman parameter 
integration turns out to be the only one needed in the calculation of general cut integrals 
in section 3.2. This is why we show it explicitly here. 


1 


— dx- 


1 


(3.14) 


(A'i,„iA“A»)(Q„4A“A») Jo (((l-a:)A- + M3).iA“A»)2 

Performing hrst the A and A integrations, we get 

1 

/ dx— ---(3.15) 

Jo {{l-x)K + xQ)^ 

The result of the x integration in the kinematic regime corresponding to the channel under 
consideration, i.e., where Kf > 0 and both iF| < 0, < 0 is then: 


^ Ta 


L A 2Q^ - ^ 2{-K2 ■ Ks - Kf) - v/a A 

^\2Q^ + VA 2{-K2 ■Ks-Kf) + VA)' 


(3.16) 
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where A denotes the discriminant that arises from the qnadratic equation in the Feynman 
parameter, given by 

A = {Klf + {Klf + {Klf - 2KfKi - 2iF|iF| - 2iF|iFf. (3.17) 


3.2. General Cut Integrals 

As discussed in section 2, calculating all one-loop Af = 1 amplitudes is equivalent 
to hnding the rational coefficients of boxes, three-mass triangles and bubbles. The box 
coefficients can be computed by quadruple cuts. In principle, the three-mass triangle coef¬ 
ficients can be computed using triple cuts and the bubble coefficients require double cuts. 
In practice we find that all bubble and three-mass triangle coefficients can be computed 
from double cuts in a simple and systematic manner. It turns out that the separation 
of the three-mass and bubble integral coefficients is easily done because of the striking 
difference in the form of their cuts; see (|3.12|) and (|3.16|) . 

In the remainder of this section, we explain how to perform general double cut integrals 


(ES). 


First, find spinor-product expressions for the two tree-level amplitudes that form the 


integrand of (|3.I|) . These are tree-level amplitudes of gluons with two fermions or two 
scalars. Two recent techniques allow the calculation of those amplitudes: MHV diagrams 
and recursion relations p8|j2^. However, at this point all we need is that they 


are rational functions in the spinor products. 

Recall that the measure in (|3.I|) is dn = + £2 — Pij), 

where Pij = pi Ppi+i +... +Pj- Using the last delta function to perform the £2 integration 
and (^3) to write the measure over £ 1 , we find 




tdt{\,dX}[X,dX]S^+\tX„\iP,f - 


(3.18) 


We denote £i by £ when there is no possibility of confusion. Recall that £aa = tXaXa- 
G{X, A, t) is the function that arises from the product of the two tree-level amplitudes in 
( p.I|) . A simple observation that helps in actual calculations is that in order to obtain 
G(A, A,t), one has to write expressions of the form (•,^ 2 ) or [•,^ 2 ] in terms of £ ~ £ 1 . A 
systematic way of doing this is by using the following identity: 


{•,h)[£2jl] {•\£2 \£i] HPrjlh] 

[£2^£i] [£2^£i] [£2^£i] 


(3.19) 
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A similar identity is valid for [•,£ 2 ]- The factors [^ 2 ,^ 1 ] and (^ 2 ,^ 1 ) all pair up in the end 
allowing for the use of the vector form of £ 2 - This happens because the product of the 
amplitudes must be invariant under the scaling zXi^ and 

Going back to the integral (|3.18|) , let us perform the t integration by using the delta 
function, 

(A, dX)[X,dX] 


a 


■ = p 


{P^^XaXaY 


■G A, A, 


P.2. 


P 'ad \ \ 

ij 


(3.20) 


Let us assume that the A dependence in the denominator of G is simpler, i.e. it has 
fewer factors than that of A. If the opposite were true, we would use the conjugate of the 
discussion that follows. 

Since A and A are independent homogeneous coordinates on two CP^ one must require t 
to transform as t > {wz)~^t when (A, A) ^ (tcA, zX) so that £aa = tXaXa remains invariant. 
For the integral (|3.20|) to make sense, it must be the case that ^(A, A, t) is invariant under 
the scaling (A, t) —> {zX,z~^t). This ensures that G in ( p.20|) has degree zero in A. This 
implies that it can be written as a sum of terms of the form 


n. {immAAiJ] 


9(V- 


(3.21) 


Each term has degree zero in A. The function g{X) contains all other terms that do not 
depend on A. 

There are two ways to proceed at this point. One is based on the introduction of several 
Feynman integration parameters. We End that this method becomes very cumbersome as 
the number of Feynman parameters increases. The second approach keeps the number of 
Feynman integrations to be at most one, and it has the advantage of leading to a clean 
separation of bubble, three-mass triangle and box coefficients. We discuss both approaches 
because they might be useful in different situations. 


Feynman Parametrizations 


We want to transform the denominator of (|3.21|) by replacing every factor 1 / [Ai , £] 
by —{Ai,£)/{£\Ai\£] and then using Feynman parameters to combine all factors in the 
denominator, including the one in (|3.20|) , into one of the form 1/{£\T{xi ,..., Xm+ 2 )\£]^~^‘^■ 
The integral to be performed has now the form 



(A,dA)[A,dA] 

{£\T{x^,...,Xm+2)\£r^^ 


Yl[AiJ]g{X) 

1=1 


(3.22) 
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where we have absorbed all into ^(A). Recall that the total degree in A of the 

integrand must be —2. Therefore, there are m factors in the numerator containing A. 

It is not difficult to hnd an analog of i.e., to write the integrand as a total 

derivative in A. Indeed, one can prove by induction that 


[A dX] e\[ri 

nLi mM 


= [M di] 




- A,-)! 


\k=0 


Here, rj is an arbitrary but hxed spinor and 
9k[xs\ = 


(m + 1 — + 1 — k) 


9k [Xi 


[rj e] 




mv] 


j + l-k 


(3.23) 


Xi^ ■ ■ ■ Xif^ 


with Xj = 


il<i2<---<'ik 


{i\T\A, 


(3.24) 


Now we can proceed as explained in section 3.1. The idea is to realize that (|3.23| ) 
is valid for all values of A except those for which there is a pole. The hnal value of the 
integral comes entirely from those poles as in (|3.9|) that led to (|3.11|) . 

The complication with this approach is that the hnal result is expressed in terms of 
several Feynman parameter integrations. In practice, if the number of Feynman parameters 
is two or more then the integration procedure is cumbersome. 


Simple Pole Expansion 

In order to avoid the proliferation of Feynman parameters we propose a second way of 
treating the integral (|3.20|) . The idea is that before doing the integral, we should separate 
the denominator factors with A as much as possible, at the cost of more terms. Where 
there is a product [a £] [b £] in the denominator, multiply both nu m erator and denominator 
by [a b]. Then apply Schouten’s identity, 

[i j] [k 1] = [i k] [j 1] + [i 1] [k j], (3.25) 

in the numerator with another factor [c 1] (which must exist by homogeneity). We then get 
two terms with [a P\ or [6 P\ in the numerator, cancelling one of the denominator factors. 
The result, in terms of A^, is a denominator of the form (£|(5^|£][H P\ in every term. 
Factors of the form {(-\Qr\(] can be treated in the same way by writing 

= [Qr^l] (3.26) 
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where Xq^ „ = -(Qr)adA^. 

Using this procedure to split poles in A and given that the integrand has degree —2 
in A, one hnds in the end only two possible kind of integrals 


Ia 


{\dX)[X^dX] ^ 

{P^tXaXaY 


Pb = 


(A,dA)[A,dA] [B4] , . 

{P^Xj^ mr\l] ^ ^ 


(3.27) 


For the second class of integrals we can apply the splitting procedure once more to 
split the product {i\P\i]{(.\Qr\(]- Then we hnd an integral of the form Ta and one of the 
form 


Pc = 


(A,dA)[A,dA] 

{l\Pm\Qr\i\ ^ ^ 


(3.28) 


Here H{X) is used to denote a generic function of i and independent of A. 

As anticipated in section 3.1, the two kind of integrals, i.e., Ta and appeared in the 
calculation of the bubble and three-mass triangle integral discontinuities. Therefore one 
can repeat the same procedure for their computation. In particular, Xa does not require 
any Feynman parameters and produces a rational function. On the other hand, Xc only 
requires one Feynman parameter and produces only logarithms. 


Canonical Decomposition 

The decomposition of the cut integral into integrals of the form Xa and Xc has a 
very useful byproduct. Since Xa produces only rational functions and Xc produces only 
logarithms, it is easy to conclude that the bubble coefficient is given by 


C2 = 


(A, dA)[A, dX] 

{Pt^XaXaY 


H{X). 


(3.29) 


Here we have used the fact that for a given unitarity cut, there is only one bubble integral 
with the corresponding branch cut. 

In the calculation of Xc one has to introduce a Feynman parameter x to write the 
factors with A in the denominator as {^\xP -t- (1 — x)Q\l] . The integration over A and A 
produces two different kind of terms. The hrst comes from the pole {l\xP -f (I — x)Q\rj\. 
The second kind comes from the poles in H{X). This is explained in detail in section 3.1 


where (|3.II|) is computed. There the two kind of terms are shown explicitly. 
Let us write the contribution to Xc from the hrst term in ( |3.II ), 

P 1 

dx -—;;- 7 ^- . „,^ H{X{x)) 


{xP + (I - x)QY 


(3.30) 
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where A(a;) is the solution of {l\xP + (1 — x)Q\rj\ = 0. 

It is easy to see that the contributions to Xq from the poles of H[X) will only have 
linear factors in x in the denominator. This implies that upon integration in x they can 
only produce logarithms of rational functions of the kinematical invariants. These simple 
functions come from the discontinuity of one-, two-, and three-mass scalar box integrals. 

On the other hand, the term given in (|3.30|) can produce a more complicated object. 
There are again two cases: if the discriminant of the quadratic equation [xP-\-{1 —x)QY = 
0 is a perfect square, then we hnd more one-, two-, and three-mass scalar box contributions. 
If the discriminant is not a perfect square then we either have the contribution of a three- 
mass triangle or a four-mass box integral. 

Recall that the coefficient of scalar box integrals can be computed very efficiently by 
using the quadruple cut technique introduced in p. This implies that we can ignore all 
those contributions and concentrate only on the three-mass triangle coefficients. 

The discriminant of {xP -1- (1 — x)Q)'^ is given by 




A = 4{{{P-Q)-Qy-{P-QyQ 


(3.31) 


Therefore, we are only interested in integrals that produce a discriminant of the form 


( CT) , i.e. 


Asm - (Kfy + {Kiy + {Kiy - 2 KIkI - 2x1x1 - ^K^X. 


3- 


(3.32) 


From the calculation of the three-mass triangle cut in section 3.1 we know that if 


P = Xi then in order to produce (|3.32|) we need Qad = ;^(-Ri)ad + (-^ 3)1 

Once we have identihed the integral that has the discriminant of a three-mass triangle 
scalar integral we have to perform one more decomposition. The idea is to expand the 
integrand in simple fractions, as a function of x, until we hnd a term of the form 


^Sm 


dx 


{xP + (1 - x)Qy 


(3.33) 
This proce- 


whose coefficient we can identify with the three-mass triangle coefficient Cg™'. 

(3*3) 

dure is applied in detail in the calculation of the cut C^s ^ in section 4.2. 

We illustrate all the features of this procedure for calculating bubble and three-mass 

triangle coefficients in the next section with the example of the six-gluon amplitnde. One 

particularly challenging technical point is that generically one might expect poles in the 

denominator of the form {i\PQ\£). This is a quadratic equation for A. We explain how to 

(3) 

deal with this in the calculation of the coefficient €^^2 in section 4.2. 
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4. Example: All Af = 1 Next-To-MHV Six-Gluon Amplitudes 


There are several pieces missing in the calculation of the six-gluon next-to-leading 
order scattering amplitude in QCD. They are: next-to-MHV MHV and next-to- 

MHV Recently important progress has been made for the hrst class of amplitudes. 

In 1^, the amplitude A(l“, 2“, 3“, 4+, 5"*“, 6 "*“) was presented. In [^, all scalar box 
coefficients of the remaining helicity conhgurations were computed using quadruple cuts. 
It is the aim of this section to present all the remaining coefficients and some new forms 
for known ones. These results complete the next-to-MHV A'^^^ piece. 


4 . 1 . First Configuration: A{1 ,2 ,3 ,4+, S"*", b"*") 


This amplitude has been computed in P^. Here we rederive it to illustrate our 


integration technique. Our result will agree with but emerge in a slightly different 
form. 

The hrst observation is that all the box and triangle coefficients vanish. This is easily 
seen by examining the helicity assignments in all possible distributions. 

Thus every double cut simply gives the coefficient of the associated bubble integral. 
The nonvanishing cuts are 034 , O'er, 0234 and 0345 . Among these, cuts 0234 and 6*345 
are mapped to each other by the permutation of indices PQ,:le->6,2<->-5,3<->'4 plus 
conjugation, while cuts 634 and Cqi are invariant under P^. There is another permutation 
: 1 <->■ 3,4 6 under which 6234 and 634 are mapped to 6345 and Cqi respectively. So 

there are only two independent integrands which are given as 


and 


634 = / dM^’^""(^i,5,6,l,2,£2)A^'^^"((-£2),3,4,(-£i)) 


= / dfi 


(3|1 + 2|6]^ 


(£i|1 + 2|6][£2 4] 


[6 1][12](3 4)(5|6 + 1|2]P|i2 (£1 5)[£2 3] 


+ dn 


(1|5 + 6|4]' 


(1|5 + 6|£2](3 El) 


(5 6)(6 1)[3 4](5|6 + l|2]P|6i [2 ^2](4 0) 


6*612 = 


dp- 


(3|P345|e] 


{«l|f345|e](3 <i) 


(4.1) 


(4.2) 


[6 1]112]{3 4)(4 5)P|^5 {<i|P345|2](5 6 ) 

Cqi 2 - When we calculate the integrand, we need to use the tree-level amplitude of four 
gluons and a pair of fermions and complex scalars. These amplitudes are given in Appendix 

B. 
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Now we will do the integration. To demonstrate our method, we will do one integration 
(for the cut Cs 4 ) in detail and then simply cite the other three results, which are obtained 
similarly. 


The Cut Cs 4 : 

The integrand is given by (|4.1|) . There are two terms which are mapped to each other 
under so that the cut C 34 is invariant. Thus we can focus on the hrst term only. In 
the hrst term, multiply numerator and denominator by {ii I 2 ) and perform the t integral 
to get 


rC) _ 
'-^34 — 


(3|Pl2|6l"P|4 


|6 11|1 21(3 4){5|P6i|2]P|: 


d\(,^)[\j dX- 


{c\Pi2\mi 3) 


612 


(£i5)(£i4) {C\P^4\cf 

(4.3) 


We can rewrite this as 



=C J (A,, dX,,)[dXt, 


3) [vh] 

V (^i5){£i4) (4|i"34|^i](^i|i"34|r7] 


(4.4) 


where we have dehned the constant 

(3|Pl2|6]V|4 

|6 1]11 2](3 4){5|P6i|2]P|i2' 

At this stage, we have three poles: \li) — |4), \li) — |5), and |£i) = |P23|h]- However, 
since \ri\ is an arbitrary spinor, we can choose it to be \rj\ = |4]. It is easy to see that after 
making the above choice we reduce the integrand into 


=c 


(A^j dXi^)[dX(,^ 


(£i|Pi2|6][4 C 


[3 4]{C 5){ii 4){h\Ps4\ii 


(4.5) 


where only one pole |£i) = |5) gives nonzero contribution (for the pole \li) — |4), since the 
factor [4 £ 1 ] appears in the numerator, the residue is zero). Reading out the residue we get 


[4 5](5|P345|61(3|P345|6]" 

(4 5>[6 1][12]{5|P345|5]{5|P345|2]P|43 


(4.6) 


Note here that the result of the integration is the residue with an extra minus sign. How¬ 
ever, the coefficient of a bubble will be just the sum of the residues at the poles. 


The Results 
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Since for the other integrations, the procedure is exactly same as above, we list our 
results directly. The coefficient in the cut Cqi 2 is given by 

(3|P612|6]" 


C2:3;6 


■ |6 1]11 2]{3 4)(4 5){5|P612|2]P|i2 V (51^1215] 
The coefficient in the cut (7234 is given by 


/(3|P5|6]P|i2 ^ (31^612^2^61216] 


(2|P612|2] 


C2:3;2 — “ 


{l|^’56l|4]• 


(5 6)(6 1>l2 3]l3 4]{5|P56i|21P|6i V (2|P56i|2] 


(4.7) 


( Wm]Pli , (1|P56iP5P56i|4] \ , , 

I (2|P56i|21 {5|P56i|5| / 


The coefficient in the cut (734 is given by 

[4 5](5|P345|6](3|P345|6]' (2 3)(1|P234|2](1|P234|4]' 

(4 5)[6 1][1 2](5|P345|5](5|P345|2]P|45 [2 3](5 6)(6 1)(2|P234|2](5|P234|2]P|34 

(4.9) 


Finally, the coefficient in the cut C^i is given by 

[5 6](5|P56i|4](1|P56i| 4]' (1 2)(3|P6i2|2](3|P6i2|6]' 

(5 6)[2 3][3 4](5|P56i|2](5|P56i|5]P5m [1 2](3 4)(4 5)(51^61215](21^61212]PI 42 

(4.10) 

It is easy to check that the sum of all these coefficients is equal to the tree-level amplitude 
as required by the divergent behavior discussed in Section 2. 


Comparison with Known Results: 


The same amplitude has been calculated in |^, where the result was given by 


with 


A —aiKo[s6i] + a2Ko[s34] 


1 

2 


, Lo[s345/s6i] , Lo[s234/s34] 

0i-h O 2 - 

^61 'S 34 


,, Lo[s234/'S6i] , , Lo[s345/s34] 

+03 -h O 4 - 


■§61 


■§34 


(4.11) 


a, =02 = 

_ (3|T345|6]^(3|P345P3452 — ^3452^34516] 

(5|P345|2][6 1][12](3 4)(4 5)P|45 
_ (3|P345|6]^(3| — 5 P 345 P 345 -t- T 3455 T 345 I 6 ] 
(5|P345|2][6 1][1 2](3 4)(4 5)P|45 

7 _ (11-P 234 I 4 ]^(11^2342^234 — -P 234 -P 2342 I 4 ] 

(5|P234|2][2 3][3 4](5 6)(6 1)P|34 

7 _ (1|-P234|4]^(1|5P234-P234 ” -^2345^23414] 

(5|P234|2][2 3][3 4](5 6)(6 1)P|34 
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(4.13) 


The functions Kq and Lq are defined by 

1 


Ko(s) = 


«)■ 


Lnfr) = 


ln(r) 


e(l-2e)^ ^ 1-r 

The f un ction Kq is proportional to the bubble integral, and the function Lq is related to 
the Feynman parameter integral for a two-mass triangle integral. These two functions are 
in fact related by the identity 

Lo[si/s 2 ] Ko[s 2 ] — Ko[si] 


S2 S2 — Si 

Therefore ( [4.111) can be brought to the following form: 
1 6 i 1 63 


+ 0 ( 6 ). 


A = \ ai 


1 

+ 2 


2 S 61 — S 345 2 S 61 — S 234 

hi , 64 


Ko[s6i] + ( O 2 


1 


1 


—h - 

■sei ~ 'S345 S34 — S345 


Ko[s345] + 


2 S 34 — S 234 2 S 34 — S 345 


(4.14) 


Ko[s34] 


- + - 

S 34 — S 234 Sgl — S 234 


(4.15) 

Each quantity in parentheses corresponds to one of the bubble coefficients in (|4.7|) - (^1.10|) . 
It is easy to check that our results agree with (|4.15|). 


4 . 2 . Second Configuration: A{1 ,2 ,3+, 4 ,5"'“,6"'“) 

This conhguration has the following Z 2 symmetry: Pq, i 7 — i plus conjugation. 
With this helicity assignment, we have box, triangle and bubble contributions. The box 
part is easy to calculate by quadruple cuts. 

The nonzero box contributions come from both two-mass hard and one-mass box 
integrals. For the two-mass hard box integrals f = 2,4, 6, the coefficients are 


2m h _ 

*^4:2:2 — 


2m h _ 

*-4:2;4 “ 


2m h _ 

*-4:2;6 — 


p2 p2 

-'61-'456 


(4|P456|3]'(4 6)[3 1] 


(6|P456|1] (4 5)(5 6)[1 2][2 3] 


p2 p2 

-'^23-' 612 


[2 6](4|P612|2](4|P612|6]' 


(4 5)[6 1][1 2](5|P612|2](3|P612|2]' 


(4.16) 


p2 p2 

-'^ 45-'^561 


(1 5)(1|P56i|3]'(5|P56i|3] 


2 [2 3](5 6)(6 1)(5|P56i|2](5|P56i|4]" 

It is easy to see that and are mapped to each other under while is 

invariant. For the one-mass box integrals i == 5, 6, the coefficients are 


Im 

C4;5 


„lm _ _ 
'-4;6 ~ 


p2 p2 

-'23-' 34 


{1|J’234|21(1|P234|31- 


(5 6){6 1)(4 2)2(5|P234|2]P|34 


p2 p2 

-'^34-'^45 


(5|P345|6](4|P345|6]^ 


(4.17) 


[6 1][1 2](3 5)2(5|P345|2]P32 


345 


19 

























which are mapped to each other under P^. 

For the triangle part, as we argued in general, we need only pay attention to the three- 
mass triangle part, For this case, there is only one function with the distribution 
(23|45|61). We can calculate the coefficient by triple cut in principle, but we choose to 
read it out by a corresponding double cut integration which we will evaluate presently. 

For the bubble part, we have following cuts: three particle channels Ci 2 ^^ Cqi 2 and 
1^234; two particle channels (723, ^' 34 ? ^*45 and (Fgi. Among them, the pairs ((7612,(7234) 
and ((723,(745) are exchanged under Z 2 symmetry while others are invariant. So in total 
we have five independent double cuts with the following integrands. 

(4|P456|3]" [3£i](£i 4) 

[1 2] [2 3](4 5)(5 6)P|56 [1 6) 

(1|P56i|3]" [3£2](^2 mil] 

[2 3] [3 4] (5 6) (6 [4 £ 2 ] (^2 5) [2 £ 1 ] 

[5 6]4(4 2)^ _ (4 ii)[ii 3] 

(2 3)[5 6][6 l]P|6i(4|i"56i|l] [2 ^i](^i|P56i|5] 

(4|P456|3]" (4|P456|^2](4 2) (£2 1^23^^456 |4) 



[2 3](4 5)(5 6)P|^6(4|A56|1] (6|P456|^2](4 3)(£2|P23|1] 


1 (1 £2 + 4 + 5 5] [3 £ 1 ] (£1 

|P23|5]((1|6|5](2 £i) + (2|3|5](1 h))^ 

' [4 5](6 1)(2 3) (£2 + 4 + 5)2 [2 h] {£i\ 

|i"23|4] {ill 

|P6l|5](6|P45i"23|^l) 


Cs4 — dn 


[3 4](4|P612|6]^ 


[3 £2](^2|P612|6] 


[6 1][1 2]Pi,Pim\P6i2\2] [ 4 £ 2](^2 5 ) 


(3 4)(1|P56i|3]^ 


{£i 4)(1|P56i|^i 


Cei — dn 


+ 


(5 6)(6 1)P|4P|6i(5|P56i|2] {il3)[2e^] 

[16](12)2[3 5]4 [6i2]{2ii) 


[3 4][4 5](2|P345|5]P|iP|45 [1 ^2 ](^i|P345|3] 
(6 1)[5 6]2(4 2)4 (£1 1)[5£2] 


(2 3)(3 4)(2|P234|5]P61^234 {il 6)(4|P234|^2] 


1 (4|£2 + 2 + 3|3] (£i| 

|P6i|3](4 h){l £i)((4|Pi2|3](1 C) - (4 l)(£i|6|3])2 

(6 1)[2 3](4 5) (£ 2 + 2+ 3)2 

(^il 

|i"6li"23|4)(^l 5)(£i 6 )(£i| 

|i"6l|2](^l| 

-P 45 3] 


{hi)|2 3|{4 5) (t2 + 2 + 3)^ (ti|iFii"23|4)(ti 5)(ti h)(ti|iFi|2|(tiM ^ 

(4.18) 

Again, when we calculate integrands we need to use the tree-level amplitude of four gluons 
and a pair of fermions or scalars given in Appendix B. 

Now we will perform the integration. Compared to the previous subsection, this 
integration is more involved. We use two examples to demonstrate our method. 


The Cut (7234 
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The integral in ( f4.18|) may be rednced to 


C234 =C' 


j{idt}lide] 


^234 |3 e] (e l)(<'|f234|3] 

imuVf e\ {1 5){«lft34|2i 


(4.19) 


after the t integration, where £ = £2 and 

„ (1|P234|31" 

12 3]13 4]{5 6>(6 1)P|34' 

Now we can see the new featnre in the above integrand: it depends on the antiholomorphic 
variable |£] as well as the holomorphic variable \t). To simplify the calcnlation, we split 
the above integrand by multiplying numerator and denominator by (^|T234|4]. Then use a 
Schouten identity to rewrite [3 ■£]{^|T234|4] = [3 4](£|P234|^] — [3|P234|^)[^ 4], we get 


^234 ^CR 


234 


dt,\i di] 


{£ 1)(£|P234|3] 


{i 5)(£|P234|2](£|P234|4] 


[3 4] 


[4 ^]{^|P234|^] 


+ 


(^|i"234|3] A 
(^1^2341^] V 


.^( 1 ) 

■^234 


_j_ /-i( 2 ) 

I *-^234 


(4.20) 

Splitting the whole integral into two pieces not only simplihes the calculation but also 
provides a nice way to separate the various contributions. As we will see shortly, the hrst 
term will produce a pure logarithmic contribution which is related to the imaginary part 
of the box integral while the second term produces a rational function which is exactly 
the coefficient of the bubble function. This same pattern will show up in every calculation 
we meet. Since we have already found the coefficient of the box integral from a quadruple 
cut, we can neglect this term and concentrate on the rational piece only. In other words, 

this splitting is the canonical way to separate box, triangle and bubble contributions. 

( 2 ) 

Now let us do the integration term by term. For <7234 we have 


/^( 2 ) _ ID 

'-^234 — 


234 


{£ d£)[e de] 


{£ 1)(£|P234|3] 


I-P 234 I 3 ] 


{£ 5)(£|P234|2](£|P234|4] (£|P234|^]^ 


= - CPis, 


{£ l)(£|P234|3]"[r7 £] 


(4.21) 


{£ 5)(£|P234|2](£|P234|4](£|P234|h](^l -^ 234 !^] 


pole 


where ( )poie means to sum the residues of all poles, which is the consequence of the 
integration f {£ d£)[£ d£]. Choosing |? 7 ] = IP 234 II), we hnd that there are three poles giving 
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nontrivial contribntions: \l) = |5), \t) = I-P 234 I 2 ] and \l) = I-P 234 I 4 ]. Snmming np these 
contribntions we finally get 


,(2) (l|i^234|3]^ ( (1|P6|5](5|P234|3]^ 

234 [2 3][3 4](5 6)(6 1)P|34 (5|P234|2](5|P234|4](5|P234|5] 

, (12)[2 3]2p|34 (14)[3 4]2p|3, \ 

[2 4](5|P234|2](2|P234|2] [4 2](5|P234|4](4|P234|4]y' 


For Cgi 


we have 


/^(l) _ 

'-^234 —'-^-'234 


234 


j{£ di)[£ d£] 


1)(^|P234|3] [3 4](£4) 

{£ 5)(£|P234|2](£|P234|4] (^b4|^](^|P234|^] 



j{£ d£)[£ d£] 


(£1)(£|P234|3][3 4](£4) 1 

{£ 5)(£|P234|2](£|P234|4] {£\P\£f^ 


(4.22) 


(4.23) 


where in the second line we have nsed the Feynman parametrization to rewrite the in¬ 
tegrand with P = 2:^4561 — P 4 . At this stage, the integration is easy to do and given 
by 


_ / (<;|P234|3]13 4](C 1) (4|P|<;] \ 

' V(^ 5)(£|P234|2](£|P234|4] PH£\P\£]) 

(5|P234|3][3 4](15) (4|P|5] [2 3][3 4](1|P234|2] (4|PP234|2) 

(5|P234|2](5|P234|4] P2(5|P|5] + (5|P234|2][2 4] P^2\P2S4PP234m 

[3 4]"(1|P234|4] (4|PP234|4) \ 

[4 2](5|P234|4] P2(4|P234PP234|4]; 

(1 5)(1|P234|3]^(5|P234|3] / ^ PI 4 ^ -Pi ^ 

[2 3](5 6)(6 1)(5|P234|2](5|P234|4]" V PIa “ Pi PIa " Pi) 


/^(l) p2 

'-"234 ~ '-"-‘234 


=CP 


234 


dz 


(1|P234|2](1|P234|3]^ ( PlPl 

(5 6)(6 1)[4 2]2(5|P234|2]P|34 U"P|34 - PI){PIa 



(4.24) 


It is easy to see that the first logarithm is the contribntion of two-mass-hard box integral 
and the second logarithm is the contribntion of the one-mass box f un ction 


The Cut C 23 

Now we consider the cnt C 23 from which we can read ont the coefficient of triangle 
fnnction. The integrand is given by three terms. For the first two terms, the calcnlations 
are similar to those above, so we jnst list the resnlts for coefficients of bnbbles (and neglect 


22 



those of boxes). They are 


^ _J1) , J2) , J3) . 
<-2:2;2 —'-'2:2;2 ' ’^2:2;2 ' *^2:2;2’ 


'- 2 : 2:2 


'- 2 : 2:2 


(2 4)2[5 6]4 (4|P56i|5] (5|P56i|3] 


[5 6][6 Ij-Plei (4|P56 i|1] (3|-P56i|5] ( 5 IP 561 -P 23 P 561 I 5 ] 


(4|P456|3]' 


[3 2](2 1)(4|P456|1]" 

[2 3](4 5)(5 6)P|56(4|P456|1] ^ [2 1] (6|P456|1] (l|i^23|l] 
(4 6)2[2 3](2|P456|6](P|56)^ 


(6|-P456|2](6|P456|1](6|P456-P23-P456|6] 

The third term is a little involved. Dehning 

(£|P23|5]((1|6|5](2 £) + (2|3|5](1£))2 


Q = 


(£|P23|4] (^|P6l|5](£|P23i"45|6) 

61-^23 + -^23-^61 ^2 


p2 

-'23 


Q = 


p2 p2 

-'^61-'^45 


p2 

-'^23 


it is easy to see that the third integration becomes (with C = [4 5 ]( 6\)(2 3 ) 

[3 £] (l|P6i|5](£|P23|^]+P|3(l|^|5] 


) 


C. 


(3) _ 
23 ~ 


-C {£ di)[i di] 


p2 

-'^23 


(^|P23|^]" [2 


{i\Q\i]Pis 


9i^) 


(4.25) 


(4.26) 


(4.27) 


after the t-integration. Now we can nse onr method to split the antiholomorphic part in 
the denominator as 


eg* =c I {e M}[e di] 

1 


1 mi3‘2]f ( 110151 + 


L(f|i’23l<]12<l {<|Q|21 
gWP-Mi (){t\P23\m 2) 


{i 3> 


231 


1 


{i ^){i\P23Q\£) 

9{i){£m] 


(4.28) 


{i\P 2 s\i]{i\Q\i] {i\Q\2]{i\P2sQ\£) 


((ii6|5](£|P23Qi^)+p|3(i 


The hrst, second and third lines of (|4.28|) are respectively the contributions from the box, 


bubble, and triangle. For the second term, it is easy to read out 


1 / ^(£)(1£)(£|P23|5][3£] \ 

[4 5](6 1) V(^3)(£|P23Q|^)(^|P23|^]2'pole 


(4.29) 


There are hve poles giving non-zero contributions: three from g{t) and two from the factor 
(^|.P 23 Q|^)- To hnd the location of the last two poles, we use the following method. Taking 
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two arbitrary external momenta |a) and \b), since the spinor is two dimensional, we can 
represent \l) as 

\i) = {\a) + x\b)) (4.30) 

with nndetermined complex variable x.i Pntting it back into {£\P 23 Q\i), we get a qnadratic 
eqnation whose solntions are 

_ -((a|P23Q|^>) + {b\P23Q\a-)) ± {a b)y/A^ 

" 2{b\P23Q\b) 

where 

A3„ = (P|i)" + (pIj)" + (PI 3 )" - 2P|iP|3 - 2P|3P|j - 2P|jP|i 

Now we can write the residue as 


(4.31) 

(4.32) 


(3) 1 9W(l<)(<|f23|5]13<] \ 

14 5]{6i)^V' 


(4.33) 


with the following hve poles: 


|£l) =|P23|4], 1 ^ 2 ) = iPeilS], 14) = |i"23i"45|6) 

14) =|a) + x+|6), 14) = |a) + x_|6), 


(4.34) 


It does not seem useful to write the expression (|4.33|) explicitly because it is rather long 
and complicated. However, its structure is very clear and easy to implement in a computer 
program. One has to evaluate the right hand side of (|4.33| ) for the hve values of \t) given 
by ( [I.34|) and add up the obtained contributions. We hnd that the (|4.33|) is the most 
convenient way to present the answer. 

It is worth remarking that although the square root shows up in above expression, the 

hnal result 02 . 2-2 rational. A similar feature is encountered in calculating the coefficient 

of a four-mass box integral from a quadruple cut [0. 

(3) 

Now we are left with the third term of C 23 which can be expressed as 


a 


(3;3) _ 
23 ~ 


-c[ + /f3(i wm) 

(4.35) 


^ In principle we should write \t) = a(|a) -|- x|6)), but one can check that the factor a drops 
out of the final expression, so we can set a = 1. 
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with P = (1 — z)P23 + zQ. Do the integration J{£ di)[i di], and we are left with 


/^( 3 ; 3 ) _ ^ 

'-^23 


=C 


dz 


[v 


gmm] 


0 L(^l^l^](^l^l^] {i\Qm\P23Q\£) 


{{im{i\P23Q\£) + Pi3{^i){i\Q\5]) 


pole 


dz 


]P\ 


0 [{i\P\i]P^ {£v){i\Qm\P 23 Q\i) 


niA ((l|6|5](f|P23QI«> + i’l3{l fW\Q\S\) 


- pole 

( 4 . 36 ) 

where at the second line we have chosen \ri\ = \P\r]) with arbitrary rj. The advantage of 
this choice is that now every pole is independent of the Feynman parameter and we can 
evalnate dz before taking residnes of poles. The integration is of the following pattern, 




dz 


{zci + C2) 


0 (ao- 2 ^ + ai^; + a 2 )( 2 : 6 i + 62 ) 


= dz 


bi{-b 2 Ci +61C2) 


1 


0 { 0 ' 2 bi — 016162 + < 3 ^ 062 ) (^61 + 62 ) 


+ 


+ 


dz 


dz 


(62C1 - 61C2) 


(22:00 + Oi) 


2 ( 026 ^ - O 16162 + Oo 6 |) (oo-2^ + Oi2 + O 2 ) 
( 20261 C 1 - 0162 C 1 - 0161 C 2 + 20062 C 2 ) 


2 (o 26 f - 016162 + oo6|) 


(oo^^ + Oi 2 + O 2 ) 


where 


( 4 . 37 ) 


0-0 —{Q — -^ 23 )^, — 2 P 23 ■ (Q — P 23 ), 0-2 — P 23 

61 ={£\{Q - P 23 )\il 62 = (^|P 23 |^] Cl = mq - P23) 


C 2 = HP'. 


( 4 . 38 ) 


23 


We have split If into three terms. Among them, the hrst two terms give the imaginary 
part of the box integral while the last term is exactly the cnt contribntion of the three-mass 
triangle fnnction. Let ns dehne the fu n ction 


Pi(oj, 6j, Cj) 


(20261C1 - 0162C1 - 0161C2 + 20062C2) 


( 4 . 39 ) 


2(026^-016162 + 006!) 

We wonld like to point ont that there is an important snbtlety with the above procednre.li 
It is correct only if the denominator 026^ — 016162 + 096! does not vanish. However, it 
can be shown that this denominator vanishes for \£) satisfying {£\P23Q\£) = 0 which is one 


® This subtlety was overlooked in previous versions of this paper. We would like to thank P. 
Mastrolia for stimulating discussions pointing to this issue and R. K. Ellis who independently 
found the problem and informed us. 
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of the poles in eq. ( [4.35|) . This means that we have to redo the simple pole expansion 
assuming that a 2 bl — aibib 2 + a^b^ = 0. In fact, we need only the term proportional to 
(apz^+aiz+aa) because Only such a term contributes to the three-mass triangle upon the 2 ;- 
integration. Redoing the simple fraction expansion we hnd that is multiplied 

by the function -R 2 («j, bj, Cj), where 


R2{(ij, bj^ Cj) 


{(l2b2Ci + a2biC2 — ni^>2C2) 
b2{2a2bi - 0162 ) 


(4.40) 


and all a’s, 6’s and c’s are as before. Thus, we hnally have the coefficient of three-mass 
triangle as 


„ 3 m _ 

Co.o.n.o — 


[4 5](6 1)(2 3) 


X 


E 

i=l 

+ 


{i k) 


g{e){e\Q\3]Ri{aj,bj,Cj) 

{ig){i\Q\2]{e\P2,Q\e) 


((l|6|5](£|P23Q|^)+i^l3(l ^)m\b]) 


l^ii 


1 


(4.41) 


[4 5] (6 1)(2 3) 


X 


E 




ff(£)(£|Q|3]R2(a„6„c,: 

{iJ^{£\Q\2]{i\P2sQ\i) 


((l|6|5](£|P23Q|^)+i^l3(l ^)m\b]) 


l^ii 


Here, in the hrst term, the summation is over the poles 


\h) =|P23|4], 1 ^ 2 ) = |P6i| 5], 14) - |i"23i"45|6), 

14) =ig|2], 14) = Ih). 

In the second term, the summation is over the remaining two poles 


(4.42) 


14) = |a)-Hx+|6), 14) = |a)-Hx_|6). 


(4.43) 


We can choose \rj) properly to reduce the number of poles further, but we do not do so 
here. It was checked numerically that eq. (|4.41|) is indeed independent of \r]). 


The Results 

Now we list the coefficients of the amplitude H(l“, 2“, d"*", 4“, 5"*", 6+). The box co¬ 
efficients are given by (|4.16|) and ( [4.17|) . The three-mass triangle coefficient is given by 
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eq. (|4.41j) . For bubbles, we have six coefficients. From cuts in three-particle channels we 
have 


(1|P234|3]" ( (1|6|5](5|P56i|3]^ 

[2 3][3 4](5 6)(6 (5|P56i|2](5|P56i|4](5|P56i|5] 

, (12)[2 3]^P|6i I (14)[3 4]^P|ei ^ 

[2 4](5|P56i|2](2|P56i| 2] [4 2](5|P56i|4](4|P56i|4]; 

(4|P612|6]" ( (2|1|6](4|P612|2]^ 

[6 1][1 2](3 4)(4 5)P|i2 (5|P612|2](3|P612|2](2|P612|2] 

I_ [5 6] (4 5)^P|j^ 2 _I_ [6 3](3 4)^Pg\2 ^ 

(3 5)(5|P612|2](5|P612|5] (3 5)(3|P612|2](3|P612|3]; 

_ (4|i^456 _ ( 16 31(6 4) _ (1|2|3](4|Pi,3|11 A 

[12][2 3]{4 5){5 6) b6|Pm|6](6|J’i23|ll {l|J’i23|l](6|Pi23|llP&3/ 

From cuts in two-particle channels, we have 


(4.44) 


[3 5](5|P612|6](4|P612|6]^ (2 4)(l|P56i|2](l|P56i|3]^ 

[6 1][1 2](3 5)(5|P612|2](5|P612|5]P|i2 (5 6)(6 1)[2 4](5|P56i|2](2|P56i|2]P|6i 

(4.45) 

and 

(2 4>"[5 61‘‘ (4|P36i|5] (51^38113] 1 

|5 6][6 {4|P36i|1] {3|P36i|5] (5|P56i -P 23 -F 561 I 5 ] 

(4|P456|3]" ( [3 2](2 1>(4|P456|1]" 

|2 3]{4 5)(5 6)P|53(4|P456|ll ( |2 l]{6|P456|l](l|P23|ll 

(4 6)^(2 3](2|P456|6](P|33)" \ 

(6|P456|2](6|P456|1] (61^456^23-^45616] / 

1 V (if f 

[4 5](6 1)^V {e^){i\P23Q\i) 

The coefficient C 2 : 2;4 can be obtained C 2 : 2;2 by a flip symmetry, while the coefficient C 2 : 2;6 
can be expressed as Atree — Xlothers'' 2 :r;i from the divergence equation (|2.(j|). We have 
obtained an analytic expression for C 2 : 2;6 using the techniques of this paper and checked 
numerically that ( |2.8|) is satisfied. The analytic expression is rather long, so we omit it 
here for brevity. 


4 . 3 . Third Configuration: A{1 ,2+,3 ,4+,5 ,6+) 

Now we move to the last helicity configuration 24(1“, 2"*“, 3“, 4"*“, 5“, 6"*"). It has the 
largest symmetry Xq generated by Pq, : z —> z -|- 1 plus conjugation. Because of this, we 
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need to calculate just one coefficient for each type of function and act on it by to obtain 


all the others. The box coefficients are 


2m h 
^4:2;1 “ 

p2 p2 

-'^345-'^56 

(l|i’345|4t 

(1|P345|5](6|P345|4] 

2 

(1 2)[3 4](6|P345|5]- 

' (2|P345|5](6|P345|3] 

and 

Prn _ 

ffi;! — 

p2 p2 

-'^45-' 56 

(5|P456|2]" 

(4 P 456 2] (6 P 456 2] 

2 

11 2112 3] (4 

(4 +456 1] (6 +456 3] 


(4.47) 


(4.48) 


For the triangle integrals, there are two two-mass triangles related by Pa- For the 
bubble integrals, the orbit of the cut in a three-particle channel contains three elements, 
while the orbit of the cut in a two-particle channel contains six elements. The representative 


integration we will perform is the following: 

(5|Pi23|2]" (5|Pi23|^2](^2|Ti23|2] [2 £2](5 4) 

[1 2][2 3](4 5)(5 6)p223 (6|Pi23|^2](^2|Ti23|l] [3 £2](4 4) 

_ [4 6]4(13)2[12] [£i2](£i3) (3 5)4(1 2)[2 6]^ (1 £2)[6 £ 2 ] 

[4 5][5 6](3|P456|6]P|56Pf2 [^1 1](^i|A 56|4] (3 4)(4 5)(3|P345|6]P|45^f2 (2 4)(5|P345|^2] 

1 (51^2 + 3 + 4|4] (5 £i)(1|Pi2|^i] (^i|Ti2|4] ((5|6|4](1 i^) + (1|2|4](5 ii))^ 

(12)[3 4](5 6) (£ 2 + 3+ 4)2 (6 £i)(2|Pi2|^i] (51^34^21^1) (£i|P56|4](£i|Pi2|3] 

(4.49) 

Now we perform the integration for these two cuts. 


The cut (7123 

After the t-integration and splitting we end up with 


<7i23 -C 


^f23(5|i^l23|3] 

(6|Pi23|3] 


{£ de)[e dh 


{^\p^2^\m 5) 


1 


I-P123I2] 


+ 


[2 3] 


(£|Pi23|1](^ 4) (£|Pi23|3] V (£|Pi23|^]" (^|A23|^][3 £] 


+ (7 


^1^23 (5 6) 

(6|Pi23|3] 


(£ dl)[l di' 


{e\Pi23\2]{e 5) 1 


(^|i"l23|2] 


+ 


(6|Pi23|2] 


(£|Pi23|1](^ 4) {i 6) {e\P^2s\ef (^I^123|^](6|Pi23|^] 

(4.50) 

n mn ot/a c\/c id -2 ■ By our Standard method, it is clear that among these 
[1 2] [2 3] (4 5>(5 6)Pj"23 

four terms, the hrst and third terms will give rational functions and the second and fourth 
terms will give the logarithmic functions involved in the box contributions. Carrying out 


with C — 


( 5 |Pi 23 | 2 ]- 


the integration for the hrst and third terms we read out the coefficient of bubble C 123 as 

(5|Fi23|21" ( (4|P„3|2f(4 5)(5 6>[4 6] 

[1 21|2 31(4 5>{5 6)P2,3 ({4|Pi33|ll(4|Pi23|3]{4|Pi23|4](4 6> 

(1 5)11 21^Pf33(5|Pi23|l] (S|Pi23|3] (5 3)12 31^Pf33 

[1 31(4|Pi23|l]{6|Pi23|l]{l|Pi23|ll {6 |Pi23|31 [3 1]{4|Pi23|31{31Pi23|3] ■ 

(5 6) (5 4)14 61(61Pi23|2f \ 

{6|P,23|31 (6 4)(61P,23|ll(61Pi23|6] j 
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We can do similar calculations for the second and fourth terms and it is easy to check that 
they produce the imaginary parts of box integrals. 


The cut Ci 2 

Here we consider the cut Ci 2 - The integrand consists of the three terms. Integration 
of the hrst two terms can be performed by the same procedure as presented earlier in the 
paper. Therefore, we will only state the results. The hrst two terms contribute to box 
coefhcients, which can be neglected for our purposes, and to the bubble coefficient C 2 : 2 ;i- 
The corresponding contributions are 


,(i) [4 6]"(1 3)^ _ (3|P456|4](4|P456|2] 

^■^4 [4 5][5 6](3|P456|6]P|5g (2|P456|4](4|P456Pi2 -P456|4] 


(4.52) 


and 


(2) _ (3 5)4[2 6]2 (1|P345|5](5|P345|6] 

(3 4){4 5>(3|P345|6]-P3t5 (5|P345|11{5|P345 Pi2P345|5]' 


(4.53) 


Now we consider the last term. After performing the t-integration, its contribution to the 
cut Ci 2 can be written as follows 


,31^ . r {i;de)iedi] ^., _,„ie2] f {i\p,^\i]{5m , (5Q1C4] \ 

J {i\Pu\ef phmit] mn)' 


where 


C = 


9{^) = 


and 


(12) [3 4] (5 6)’ 

(<5) (<|Pi3|4] ((5|8|4](1Q + (1|2|4](5Q)^ 

{( 6 ) (<|Pi2P34|5> {(’|P66|41(<|Pi2|3] 

Q=-^(PlP^-i + PhPx)- 

W 2 


Now we use our method to split the integrand in (4.54). We hnd that 


=c {i di)[e di] 


1 


gm 2 ] 


[{i\p^ 2 mli] {i\Q\i] 


(5|6|4] - 


(1 2) [4 1](5 £) 
{i 2) 


9{e)p?2{5 ^)(^|Ti2|4](£ 1 ) 


12 


{£ 2){£\Pi2Q\i) 
9{i){£\Q\2] 


{£\Pi2\mQ\i] {£\Q\mPi2Q\£) 


{{5\6m\Pi2Q\i]+p!2{5£)mm 


(4.54) 


(4.55) 


(4.56) 


(4.57) 


(4.58) 
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Among these three parts, the hrst one contributes to one of the box coefficients and can be 
neglected for our purposes. The second term gives a contribution to the bubble coefficient 
C 2 : 2 ;i- Performing the ^-integration, we find this contribution to be of the following form 






9((’){5<){flPi2|4)(2<] 


(3 41(5 6) ^ V' ' {I m\PimW\Pi2\l\) 

Here the locations of the six poles is as follows: 

|£i) = |P56|4], 14) = |Pi2|3], 14) - |6), 14) = IP 12 P 34 I 5 ), 


(4.59) 


14) = |a)a:_|6), |4) = |a)-1-a:+|6). 


with 


and 


x± = 


-{{a\Pi 2 Q\h) + ( 6 |Pi 2 Q|a)) ± (a b)y/Asm 

2 {b\Pi2Q\b) 


^3m = {{P!2f + (Puf + (P^ef - ‘^PhPk - ‘^P^Pie - ^P56P?2) 


(4.60) 


(4.61) 


(4.62) 


Note that the pole |£) = |2) does not contribute because of the factor [2 £] in the numerator 
of (14391) . 

Now we move on to the last term in ( [4.58|) . It will produce the contribution to the 
three-mass triangle coefficient. By introducing the Feynman parameter z, we can write it 
as follows: 


=C dz 


where 


{idi)[idi] -g{l){l\Q\2] 


((5|6|4](«|Pi2Q|«) + Pf2{5 

(4.63) 

(4.64) 


P — (1 — z)Pi 2 + zQ. 


Performing the ^-integration, we obtain 


(4.65) 


pole 


where \r]] is an arbitrary auxiliary spinor. Let us write 4] as 


\rj] = \P\7j) 


(4.66) 
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for some spinor |? 7 ). Then ([4.65|) becomes 


cjp' dz 




{e\p\e]p^ {en){e\Q\i]{e\Pi2Q\e} 


({hmmPiMt)+ piz{^ 


- pole 

(4.67) 


The Feynman integral of this type was considered before. It produces the functions 
Ri{aj,bj,Cj) and R 2 {aj,bj,Cj). As a result, the contribution to the three-mass triangle 
coefficient is given by 


3m 

*-3:2:2;l 


1 


E 


{Ur) 


gii){£\Q\ 2 ]R^iaj,b,,Cj) 


+ 


(12)[3 4](5 6 )^L' " {i^{£\Qm\Pi2Q\£) 

1 


((5|6|4](£|Pi2Q|^)+Pi"2(5 ^)(^|Q|4]) 


E 




(12)[3 4](5 6)^L' " {i^{£\Qm\Pi2Q\i) 

where, in this example, the coefficients a^, bj and Cj are 

do = {Q ~ Pl2)'^, Cli = 2 Pi 2 ■ {Q — P12), d2 = ^*^25 

61 = {£\Q - P,2\i], b2 = {i\Pl2\£], Cl = mQ - Pi2m, C2 - {v\Pl2\ 


In the hrst term in eq. ( 4.68 ), the summation is over the poles 


|£l) = |P 56 | 4 ], 14 ) = |Pi 2 | 3 ], 14 ) = |6), 14 ) = 1 ^ 12 ^ 3415 ), 

14 ) = 1 ^), 14 ) = IQ|1]. 

In the second term, the summation is over the poles 


(4.68) 


(4.69) 


(4.70) 


Itz) = \a} + X-\b), 14 ) = |o) + i+| 6 ). 


(4.71) 


Thus, the bubble coefficient C 2 : 2 ;i is represented by the sum 


„ _ 41 ) I „(2) , „(3) 

o2:2;l — *^2:2;1 ' *^2:2;1 ' *^2:2;1’ 


(4.72) 


where *^ 2 ^ 2 ;! c^^.i are given by ( [4.52|) , ( [4.53|) and ([I.59|) respectively. The three- 

mass triangle coefficient C 3 l^. 2 .i is given by ( |4.68|) . 
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4.4- The amplitude A{1 ,2 ,3 , 4+,..., n"*") 


The next-to-MHV n-point amplitude with all negative-helicity gluons appearing con¬ 
secutively has been computed in [^. We perform the calculation here as well, as an 
illustration of our procedure. The solution will emerge in a different form. 

From the viewpoint of our discussion in Section 2, there are no box and three-mass 
triangle contributions, hence no one-mass and two-mass triangle contributions. All we 
have are bubble contributions which can be calculated by double cuts. There are two 
nonvanishing double cuts, which we denote by I and II. 


I 

II 


dfi AL{£i,k + 1, ...,n, 1, 2, £ 2 )Aji{£ 2 , 3,4,..., k,£i) 4 < k < n - 1 

dn Al(£i, 2,3, ...,k,£ 2 )AR{£ 2 ,k + 1,..., l,^i) 4<k<n-l 


(4.73) 


However, these two kinds of cuts are mapped to each other under the permutation Pq, : 
{n — i + 4) ^ i and the identihcation k = n — k + ?,. 

Now we focus on the case I. To do this we need to know the tree-level amplitudes 
with two fermions or complex scalars, which we list in Appendix B. First we have 


Az.(£+,(/c + l) + ,...,n+,l-,2-,4) 


+ 1-0- 0 -^ 


i—k — 2 


- E 


{n-l-j n - J>(l|ifF7+„+TK+ 


n-j 


j=0 


(h + !){* + 1 fc + 2)... (n - j|A'W|21(n - 1 - 




X 


n-j 


n-j 




(1i fc+l>{l|A''" 




fc+1 


1 ^ 2 )' 




1 


X 


^[n-t-2 — (/i:-t-l)],[n-t-3 — (fc-t-1)] 
/l|j^[’^+2 —(fc+1)] j^["-+3 —(fc+l)] I ■ ' “ 


[n+2-(fc+l)] 


| 2 ] 


1 ^ 1 )' 


{h\K^ 


[n+2-{k+l)] 

k+1 


| 2 ] 


/ I j^["-+2 —(fc+l)] j^["-+3 —(fc-l-l)] 


1 ^ 2 ), 


(4.74) 

where the last term is the special case with j = n — k — 1. In this formula, a = 0,1, 2 for 
£i ,£2 to be gluons, fermions and complex scalars. Using (^4.74|) and MHV-amplitude of A^ 
we get the integrand of cut I as 

{k k + l){12){2 3) (£i 3) 


1 = - 


d/i- 


X 


n — k — l 

E 

j=o 


n 


nr=i(**+i) {£ik+i){k£i) 

1 - j n- 


(4.75) 


-i+'k7'i2i{» -1 - iiA'k7'i2i 


n-j 


b+2] I 

n-j 
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Notice that in the above integration, only the holomorphic part |£i) appears in the inte¬ 
grand, so the integration is very easy to do, similar to the example A{1~, 2“, 3“, 4+, 5"*“, 6"*“). 
We separate the integral into the cases where j ^ n — k — 1 and the special case of 
j = n — k — 1. For the cases where j ^ n — k — 1 there are three potential poles: 
|£i) = |£i) = \k + 1) and |£i) = \k). By choosing \r]] = we get rid of 

one of them. Similarly, for the case j = n — k — 1, we take \r]] = ^^|3) to get rid of 

one pole and leave only two poles: \li) = \k) and \li) = |iF2^ ^^|2]. Adding all these pieces 
together, we hnd that the coefficient from cnt I is 


C2:fc-2;3 


(1 2)(2 3) 

nr=i(^^+i) 


n—k — 2 

E 

j=0 


{n-l-j n - j>(l|A'17’W?|3)^ 
&&{n - jii&m - 1 - ji A-yyi 12] 


X 


+ 


X 


imtM 


( 3 |Af + ii(t+ "b]{*iAENdA'ii) 

(fc+iiAy^^iii+i) 

(1 2)(2 3) {k fc+ l)(l|A|'‘~‘^A'f~^'|3)^ 

n”=i(* *+ 1 ) l|A' 2 ‘’"''| 21 {fc|Ai'‘"‘'| 2 ] 


{3|Af ^'\k]{k\K 


[k 2] j^[k 1] 
O XV o 


|1) (2|A:f ^^|3)(l|A:f ^'|2]' 


{fc|A-f-"l|fc] 

By symmetry we read ont the coefficient from cnt II as 


(214''' "'|21 


(4.76) 


C2-.k-l-2 — — 


X 


2»2 3) „ (j+5j+4)(3|K3"' '|1)- 

nfei{” + l) 4-'+"'4-'+"l{2 + 4|Af'""l|21{i + 5|A4"'|2] 

/ (llAbthi" - fc + 31(» - *= + -""'IS) 

V (n-fc + SIA'lyffiln-fc + S] 

(llAjtthln -k + 41(n -k + 4|A4"lAf-""llS) 

'~-fc + 4|Afrd4l"-* + 41 

{n-k + i n-k + aXSIAfrffiAfrd^ll) 


(n 




(1 2)(2 3) _ \n-K-r^ n -/v--e 

r=i{i i + 1) 4Al44'‘-"d4(" - + 3|Ai'rd4l21(n - k + 4\Kt,+i 

A“2] T^lk—l] |q\ 
n-k+4:-^n-k+4A/ 


[fc-1] 


X 


( -k+i]{n-k+ 41^4714^4714 
V (n-fc + 4|Af-4jfc-fc + 4] 

(2|4tyi4AlEi4|l)(3|A4td4|2] ' 

(ai/WAiaJ 


12) 


(4.77) 
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In [^1, the amplitude was decomposed in terms of functions Kq and Lq, given here 
in (|4.13|) , where the function Kq is proportional to the bubble integral and the function 
Lq is related to the Feynman parameter integral for a two-mass triangle integral. By the 
identity (|4.14|) , it is possible to convert their expression to an expansion in bubble integrals 
only, as we have done here. 


5. Summary Of Results For Next-To-MHV Six-Gluon Amplitudes 


In this section we collect our results for the next-to-MHV six-gluon amplitudes for the 
reader’s convenience. 


We dehne the following functions: 


Ri (gj , bj , Cj ) 
R2 5 bj , Cj ) 


( 2 a 26 iCi - ai 62 Ci - ai6iC2 + 2aob2C2) 
2 (g 2 ^i ~ G1&1&2 + 

{a2b2Ci -F a2biC2 - aib2C2) 

62(20261 - 0162) 


(5.1) 


5.1. A(l-,2-,3-,4+,5+,6+) 


The amplitude is given by 


A{1 ,2 ,3 , 4+, 5"^, 6"^) (c2:3;6-f2:3;6 + C2:3;2f2:3;2 + C2:2;3-f2:2;3 + C2:2;6-f2:2;6) 

(47rJ^ e 

(5.2) 


with rr = 


^_ (3|P612|6]^ _ / (3|P5|6]P|i2 , (31^612^2^61216] \ 

[6 11(121{3 4){4 5)(5|P6i2|2]Pey V (5|-P612|51 (2|P6 i2 |2] ) 

^ _ (1|P56i| 4]" _/ (l|P2|4]P|6i ( 1 | -P 561 -P 5 -P 561 14] \ 

“ (5 6)(6 1)[2 3][3 4](5|P56i|2]P|6i V (2|P56i|2] W^] ) 

[4 5](5|P345|6](3|P345|6]' (2 3)(1|P234|2](1|P234|4]' 

(4 5)[6 1][1 2](5|P345|5](5|P345|2]P|45 [2 3](5 6)(6 1)(2|P234|2](5|P234|2]P|34 


[5 6](5|P56i|4](1|P56i| 4]' (1 2)(3|P6i2|2](3|P6i2|6]' 

(5 6)[2 3][3 4](5|P56i|2](5|P56i|5]P|6i [1 2](3 4)(4 5)(51^61215](21^61212]P|i2 

(5.3) 


34 























5 . 2 . A(l-,2-,3+,4-,5+,6+) 


The amplitude is given by 


A(l“,2“,3+,4“,5+,6+) 


^r(h ) 


2Ne 


(47r'|2-e 

.Jim 


) ( 2m h r2m h , „2m h j2m h , „‘2.m h j2m h i Im rim i Im rim 
“ V'^4:2;2 -'4F:2;2 ' *-4:2;4 -'4F:2;4 “r <.'4.:2\G -'4F:2;6 ' *-4;5-'4F;5 ' *-4;6-'4F;6 


+ ^3^:2-2^i^:2-2 + C2:3;2-^2:3;2 + C2:3;6-^2:3;6 + C2:3;l-^2:3;1 
+C2:2;2-^2:2;2 + C2:2;3-^2:2;3 + C2:2;4-^2:2;4 + C2:2;6-^2:2;6) 


(5.4) 


with 


2m h _ 

*-4:2;2 “ 


2m h 
^4:2:4 


2m h _ 

*-4:2:6 ~ 


p2 p2 

-*61-* 456 


(4|P456|3]'(4 6)[3 1] 


(6|P456|1] (4 5)(5 6)[1 2][2 3] 


p2 p2 

-*23-* 612 


[2 6](4|P612|2](4|P612|6]' 


(4 5)[6 1][1 2](5|P612|2](3|P612|2]^ 


p2 p2 

-* 45-* 561 


(1 5)(1|P56i|3]'(5|P56i|3] 


[2 3](5 6)(6 1)(5|P56i|2](5|P56i|4]^ 


Im 

*^4;5 


Im 

*-4;6 


p2 p2 

-*23-* 34 


(1|-P234|2](1|T234|3]" 


(5 6){6 1>[4 2]2{5|P234|21P|34 


p2 p2 

-*^34-^^45 


(5|P345|6](4|P345|6]^ 


[6 1][1 2](3 5)2(5|P345|2]P3" 


345 


(5.5) 


„ 3 m _ 


•X 


[4 5](6 1)(2 3) 

g{l){l\Q\?,]Ri{aj,hj,Cj 


E 

j=l,2,3,6,7 




+ 


(£77)(£|Q|2](£|P23Q|^) 


■ X 


{mmP23Q\t)+ ^)(^ iq | 5 ]) 




[4 5] (6 1)(2 3) 

g{l){l\Q\?>]R2{aj,hj,Cj) 


E 

1=4,5 *- 




{lg){l\Q\2W\P23Q\l) 


((1|6|5](£|P23Q|^)+P|3(1 ^)(^IQ|5]) 




(5.6) 


35 



(l|i"234|3]^ f (1|6|5](5|P56i|3]" 

[2 3][3 4](5 6)(6 l^(5|P56i|2](5|P56i|4](5|P56i|5] 

(1 2)12 (14)[3 4]^P|6i \ 

[2 4](5|P56i|2](2|P56i| 2] [4 2](5|P56i|4](4|P56i|4]y' 

(4|P612|6]^ f (2|1|6](4|P612|2]" 

[6 1][1 2](3 4)(4 5)P|i2 1^(5|P612|2](3|P612|2](2|P612|2] 


, [5 6](4 5)^P|,, [6 3](3 4)^P|,, \ 

(3 5)(5|P612|2](5|P612|5] (3 5)(3|P6i2|2](3|P6i2|3] J 

(4|P456|3]" r [6 3](6 4) (l|2|3](4|Pi23|l] 1 

[1 2][2 3](4 5)(5 6) [(6|Pi23|6](6|Pi23|l] (l|Pl23|l](6|Pl23|l]Pf23. 


(5.7) 


(2 4)2[5 6]4 (4|P56i|5] (5|P56i|3] 1 

[5 6][6 l]P|6i (4|P56 i|1] (3|P56i|5] (5|P561 -P 23 -P 561 I 5 ] 

(4|P456|3]^ f [3 2](2 1)(4|P456|1]" 

[2 3](4 5)(5 6)P|^6(4|i^456|l] [2 l](6|P456|l](l|i^23|l] 

(4 6p[2 3](2|P456|6](P|56)^ ^ 

(6|P456|2](6|P456|1](6|P456P23-P456|6] / 

1 V L / . ^(^)(1^)(^|P23|5][3£] - 

[4 5](6 1) ^){i\P2sQ\i){i\P2s\i]. 

[3 5](5|P612|6](4|P612|6]' (2 4)(l|P56i|2](l|P56i|3]' 

[6 1][1 2](3 5)(5|P612|2](5|P612|5]P|i 2 (5 6)(6 1)[2 4](5|P56i|2](2|P56i|2]P|6i 

C2:2;4 =-Pa(c2:2;2) 

C2:2;6 =^tree ” C2:2;2 “ C2:2;3 ~ C2:2;4 “ C2:3;2 “ C2:3;6 ~ C2:3;l 

(5.8) 

The symmetric action is Pq, : i ^ 1 — i plus conjugation, i.e, ( )<->■[ ]. For (|5.6|) and 
( ^78|) , we have dehned 


<^o — (Q ~ -^23)^5 «! — 2P23 ■ {Q — 7^23)5 «2 — -P23 

61 =(£|(g-P23)|7], 62 = (7|P23|7] Cl = (i7|(g-P23)|7], C 2 = (^ 1 ^ 2317 ] 

|7i) =|P23|4], 1^2) = |P6 i|5], 14) = |P237"45|6) 

14) =|a)+x+|6), 14) = |a)+x_|6), |4) = |(5|2], |4) = lij) 


(5.9) 


(5.10) 


-((a|7^23QI4 + (^|7^23Q|a)) ± (a 
2(6|P23QI4 


(5.11) 


36 



























where 


A3„ = {Pilf + {Pif + (PI3)" - iPlPh - "iPlPl - 2P«P|i 

^ (<|P23|S1 ((1|6|S1(2 Q + (2|3|5](1Q)^ 

® ’ (<|P23|4] (<|P6l|8](f|P23P45|6) (5'12) 

Q = ■;^(^’Ii^’ 23 + PisPei) 

-^23 


Here |a), |6) and \r]) are arbitrary. 


5.3. H(l-,2+,3-,4+,5-,6+) 


This helicity conhguration has the largest symmetry, namely a Zg generated by : 


i ^ i 1 pins conjngation, so we have gronped everything into orbits. The amplitnde is 


given by 


24(1-,2+,3-,4+,5-,6+) 


rrin 


2^e 


E 2m h r2m h , \ ^ „1?ti rim 1 \ ^ „3m rSm 
C4:2;i -'4F:2;i , ^3:2-.2-.P3:2-.2-.i 


3 j=l 


i=l 


3 6 

+ C2:3;i/2:3;i + C2:2;i-f2:2; 


i=l 


i=l 


i=l 


(5,13) 
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with 


p2 p2 

^2m h _-'^345-'^56 


^4:2;! 


^4:1 — 


(1|P345|4]' (1|P345|5](6|P345|4] 


p2 p2 

-'45-' 56 


(1 2)[3 4](6|P345|5]^ (2|P345|5](6|P345|3] 

(5|-P456|2]^ (4|P456|2](6|P456|2] 


[1 2][2 3](4 6)2p2^6 (4|P456|1](6|P456|3] 


3m 

Co. 


(12)[3 4](5 6) 


X 


2,3 

1 


m 


g{l)m\ 2 ]R^{a,,h,,ci, 


{mmPi2Q\i)+ pu^ ^)mm 




E 


C2:3;l = 


(1 2>|3 41(5 6) , 

(6|Pi23|2]" 




(4|Pi23|2]'(4 5)(5 6)[4 6] 


[1 2][2 3](4 5)(5 6)Pi^ 23 \ (4|Pi23|l](4|Pi23|3](4|Pi23|4](4 6) 


(1 5)[1 2]^ 

-^123 (5|Pi23|1] 


+ 


(5|Pi23|3] (5 3)[2 3]^Pi43 


[1 3](4|Pi23|l](6|Pl23|l](l|A23|l] (6 |Pi 23|3] [3 1] (4|Pi23|3] (3|Pi23|3] 

(5 6) (5 4)[4 6](6|Pi23|2]" 


C2:2;l 


(6|Pi23|3] (6 4)(6|Pi23|1](6|Pi23|6] 

[4 6]"(1 3)" _ (3|P456|4](4|P456|2] 

[4 5][5 6](3|P456|6]P|5g (2|P456|4](4|P456 Pi2 -P456|4] 
(3 5)4[2 6]2 (1|P345|5](5|P345|6] 


+ 


(3 4)(4 5)(3|P345|6]P3i5 (5|P345|l](5|P345i"l2P345|5] 


1 


lim 




(7(£)(5 £)(£|Pi2|4][2£] 


[3 4](5 6) L 2)(£|Pi2Q|£)(£|Pi2|^] 

(5.14) 

and the other coefficients can be obtained by the symmetric action of Pq,. Here one needs 


|£l) = |P56|4], \^ 2 ) = |Pi2|3], 14) = |6), 14) = IP 12 P 34 I 5 ), 

14) = |a)+rr_|6), |4) = |a) + rr+|6), |4) = [i)), |4) = |g|l]. 


(5.15) 


O'O — (Q ~ Pl 2 )^, CLi — 2 Pi 2 ■ (Q — P 12 ), «2 — Pi 2 j 

h = {i\Q-PlM h = {i\Pl2\il Cl = (i7|(Q-Pi2)|£], C2 = (4^21^] 
a(f) = (^5) (£|Pi2|4] ((5|6|4](14 + (1|2|4](5 4)^ 

^ {£ 6) (£|Pi2P34|5) (£|P56|4](£|Pi2|3] 

Q = 4r(PkPi2 + P?2P!,6). 

^12 


(5.16) 


(5.17) 
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where 


-((a|Pi2Q|6) + (6|Pi2Q|a)) ± (a b)^/A^ 

2{b\Pi2Q\b) 

^3m = {P!2? + (Pm? + {Pi^f - ^P!2Pm - ‘^P^Pm - ‘^Pi^Pl2 

and I a), |6), \rf) can be chosen arbitrarily. 
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(5.18) 

(5.19) 


Appendix A. Scalar Integral Functions 


In this appendix we list the explicit results for the scalar integrals derived in the 
reduction procedure of 0. The expressions here are taken from [P|,|i7[|. 

The dimensional regularization parameter is e = (4 —D)/2. The constant rr is dehned 


by 


r(i + 6)r^(i-6) 
r(l-2e) 


(A.l) 


i+r+r’ 




Fig. 3 : Scalar bubble and triangle integrals, (a) One-mass triangle /a™- (b) 

Two-mass triangle (c) Three-mass triangle (d) Bubble l 2 -.r-,i- 
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i i-1 


(f) 

Fig. 4: Scalar box integrals, (a) The outgoing external momenta at each of the 
vertices are Ki, K 2 , K 3 , K 4 , defined to correspond to sums of the momenta of gluons 
in the exact orientation shown, (b) One-mass /i™- (c) Two-mass “easy” 

(d) Two-mass “hard” (e) Three-mass (f) Four-mass 
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Scalar bubble integrals: 


, ,2-e f 1 

" y (27r)4-2^p2(p_^)2 

h-.r-i = ’’r Q - ln(-4"'^) + 2^ + 0{e) 

Scalar triangle integrals: 


Is = t(47r)^ " 


‘^'^p 1 

(27r)4-2^ p‘^{p - Kif{p + KsY 


(A.2) 

(A.3) 


(A.4) 



6\r\i 


rr (Ab--(-A7~’')- 

(-7')-(-A7“b 


rSm 

^3:r:r';i 


i 



f 1 ^ 

V 


Li2 


V 1+A 




We have defined the following: 


A3 = -{Klf - {Klf - {Klf + 2KIkI + 2KIkI + 2KIkI 


51 

52 
5s 


[r] _ [r'] _ [n-r-r'] 

^i-\-r ‘'i+r+r' 

_j.W , j-h'] _ ,[n-r-r'] 

_+M _ +b'] I ,[n-r-r'] 
'^i+r I ^i+r+r' 


Scalar box integrals: 


(A.5) 


(A.6) 


7-4m 

-^4 


—t(47r)^ 


d^e 1 

(27r)4-2^ (£2 + ie)((£ _ X^)2 + 


K2)2 + A)((£ + K4)2 + A)' 
(A.7) 


We separate divergent terms from the finite pieces of interest, as explained more fully 
in section 2. We refer to the finite pieces I 4 F as finite box integral functions. 
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-M%r ? + (-tii)-' - 

iy- ot ■ ry 




M .[2] 


Lio I 1 — 


+ Lio I 1 — 




(A.8) 


r2m e _ _e _l 'l“'= _ r2™- e 

24:r;i _^[r+l]_^[r+l] _ _^[r+2]_^[r] ^2 ^ A ) \ ) I ^i-1 ) J ■‘■4F:r-i 

2rr L ("i.jLULi A - iA'l + Li A 

[r+1] [r+1] _ ,[r+2] [r] 2 I 1 + 2 I 1 + 2 I 1 

’'i —1 *^1 —1 L \ ^i—1 / \ / \ ^i — 1 / 


+ Li2 1-^ -Li. 1- 


,[r] ,[r-h2] 
''i ''i — 1 
[r-M] [r-M] 
''i-1 ''i 


1 . 9 / t 


+ - In 


2 0-1 


+ W 


(A.9) 


2m 2^r ^ rir-/[2] 1-^ , 

[2] ,[r-|-l] g2 2^ ‘' 1 - 2 ) '^ \ ''i-l 1 ^ 

o^9_ 1 L 


-*-/ j.[^+2]\—e I T2m h 

-[ C-_2 ; "T-^4F:r;i 


TZm h _ 

^4F:r:i “ 


ipi.AV L 2 V tr' 




/ Jr] \ / Jr+2]\ 

+LiJl-^jLLiJl-|^j+OW 


(A.IO) 


j6m _ 

■^4:r:r';2 


i[r+l] .[r+r'l _ . [r] , [r-|-r' + l] 

^i-1 4j rj_]_ 


+[’'+4]'\-e I ^[r+r']\—e 1/ i[r]\—e 1/ i[r-|-r'-|-l] x-e ^3m 

’ ~ 2 ~ 2 + -'4F:r:r';i- 


__llnfS- 

4F.r.r ,i [^_|_ 1 ] [^_|_^/] ^ rr_|_r/_|_l] 2 1 ,[r] 

0-1 O ^i-1 L \ 


,[r+r'] / 2 

''i+r / 




1 / 

I ^ ) 2 ’'i-l 

+ IT In 


7T + Li2 1 


,[r+l] 

0-1 


+ Li2 1 


, [r+r' + l] 
^i-1 


2 v4:iy 

/ Jr] Jr+r'-hl] \ 


fV I -Jj.1' 

''i-l ''i 


(A.ll) 


The dilogarithm function is dehned by Li2(a:) = — jo ln(l - z)dz/z. 
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Xn 


j . , T , t34.-ie 

Ll2 1 + -- -X 

ti3 - ^e 


^34 - ie\ , A , t34 - ie 
j I V \ -Xk, -, -- I In I 1 + r- -X, 


- Li2 ( 1 + ] -V \ -Xk 


ti3 - ie 

^24 ~ 

’ ti 2 - ie 


ti 3 - ie 


In ( 1 + ———Xj 
ti2 - ^e 


+ In(-Xj)(ln(ti 2 - ie) + ln(ti 3 - *e) - ln(ti4 - ie) - ln(t23 - *e))). 


(A.12) 


Here we have defined tmi = —{Km + Km+i + ... + A/-: 


r]{x, y) = 27ri['d(—Inix)'d(—Iniy)i?(Ini(a;y)) — ^?(Inlx)^?(Imy)^?(—Im(xy))], (A.13) 

and xi and X 2 are the roots of a qnadratic polynomial: 

ax^ + 6x + c + ied = a{x — xi){x — X 2 ), (A.14) 


with 

a = ^24^34, 

b = tl3^24 + ^12^34 — ^14^23: 
C = ti2ti3, 

d = t23- 


(A.15) 


Appendix B. Tree Amplitudes with Fermions and Scalars 

Here we summarize some results that are useful for our calculations.0 Take a = 2 
for a scalar and a = 1 for a fermion. (Taking a = 0 reproduces the results for all-gluon 
amplitudes, but these are not needed in this paper.) 


^ 5 3^^g) — /i2|3456 +-^1234|56 

^_ (3|1 + 2|6]^ _ f (4|1 + 2|6] \ 

[6 1][12](3 4)(4 5)P|45(5|6 + 1|2] V (3|l + 2|6]y' 
(1|5 + 6|4]^ /(l|5 + 6|3]y 

[2 3][3 4]{5 6>(6 l)P|ei(5|6 + 1|2] V(1|5 + 6|4]J 


7 


Some of these results have appeared in 



(B.l) 
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Similarly we have 


A(4-/5,5+,6+,1-,2- 

’ ^F/s) ~ -^12|3456 + +-^1234|56 

(4|1 + 2|6]^ /(3|l + 2|6]y 

^ 16 1111 21(3 4>{4 5)P|^5{5|6+1|21{3|1 + 2|6] V{4|1 + 2|6]J (B.2) 

_ (116 + 613]“ _ / (1|5 + 614] \° 

|2 31|3 41{5 6){6 1>P|6i(5|6+l|21{l|5 + 6|4] V {115 + 6|3]1 


Notice that for the split /i23|456 with a fermionic line between the two vertices, the 
amplitnde is zero (nnlike the internal glnon case). 


— /i2|3456 +-^612|345 +-^5612|34 

^_ [1 3]4(5 6)4 _ / (5|Pi23|3] y 

[1 2][2 3](4 5)(5 6)p223(4|Pi23|l](6|Pi23|3] V [1 3](5 6) J 

(6 2)4[3 4]4 / (2 1)\“ (B.3) 

^ (6 1)(1 2)[3 4][4 5](2|P612|5](6|P612|3]P|i2 V (^7 

, _ (2|P234|1]^ _ ( (2|P234|6] y 

(2 3)(3 4)[5 6][6 1]P|34(4|P234|1](2|P234|5] V (21^23411]/ 


^(^F/S’2 13"*“, 4"*“, 5 , — /23|4561 +-^123|456 +-^6123|45 

[3 4]4(5 1)4 /(5 6)y 

■ [2 3][3 4](5 6)(6 1)P|34(1|P234|4](5|P234|2] [WT)) 

_ (1 2)4[4 6]^ _ / (2|P456|4] \ “ (B.4) 

^ (1 2)(2 3)[4 5][5 6](3|P456|6](1|P456|4]P|56 V (2 1)[6 4]) 

, _ (51^34516]" _ ( (5|P345|l] r 

(3 4)(4 5)[6 1][1 2]P|45(5|P345|2](3|P345|6] V (5|P345|6]y' 


^(^F/S’2 i3"'“,4 , 5"*", — /i2|3456 +-^612|345 +-^5612|34 

^_ (4|Pi23|3]" _ f [l 3](4 6) y 

[1 2][2 3](4 5)(5 6)Pf23(4|Pi23|l](6|Pi23|3] ^41^^12313]^ 

(12)4[3 5]^ /(2 6)y (^•^) 

(6 1){1 2)13 4114 5]{21P612|51{6|P612|3]P|i2 

[5 6]4(4 2)4 / [5 1]y 

(2 3){3 4)15 6116 l\Pl^S\PMm\Pi3t\5\ V M/ 
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^(1^/5’2 ,3+,4 ,5+, ®F/s) ~ A2|3456 + -^612|345 + -^5612|34 

^ _ [1 3]4(4 6)" _ / (4|Pi23|3] y 

[1 2][2 3](4 5)(5 6)p223(4|Pi23|l](6|Pi23|3] V [1 3](4 6) J 
(6 2)4[3 5]^ / (2 1)y 

^ (6 1>{1 2)13 4114 51(2|P6i2|5](6|P8i2|31P|i2 

[5 1]4(4 2)4 /[sejV 

^ (2 3>(3 4)15 6116 Im; 

The NMHV tree-level amplitude with adjacent negative helicities is given by 

A(4+/5,5+,...,n+,l-,2-,3-/g) 

nr=3(* * + 1) -ji4"“'“"^|2](n-j - ii4"“'“"^|2] (b.t) 
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